Properties and evolution of a protoneutron star in the enlarged SU(3) model. 
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Protoneutron stars are hot and lepton rich objects formed as a result of type II supernovae 
explosion. This paper describes results of analysis of protoneutron star models constructed 
under the assumptions that strange particles are present in the core. All calculations have 
been performed for the neutrino opaque matter when the entropy per baryon is of the order 
2. The equation of state which is crucial for the structure and composition of a star has been 
obtained in the framework of the chiral SU(3) model and compared with the one involving 
derivative coupling of baryons to mesons. Of special interest in this paper is the comparison 
of the two protoneutron star models with those formed when neutrinos leak out of the 
system, described as cool, deleptonized neutron stars. It has been shown that the maximum 
mass of the hot protoneutron star with trapped neutrinos is close to 1.9M Q in case of the 
SU(3) model. However, deleptonization reduces the value of stable protoneutron star mass 
significantly to ~ 1.5M Q . 



sectionlntroduction The evolution of a nascent neutron star has to be described under the as- 
sumption of neutrino trapped matter. Numerical solutions with the lepton number Yl ~ 0.4 are 



accepted. A protoneutron star is formed 



as a result of a supernova explosion which 



leaves a hot, neutrino opaque core surrounded by a colder, neutrino transparent outer envelope. 
A protoneutron star structure and composition depend strongly on the chosen form of the equation 
of state which in turn is connected with the nature of strong interactions. However, the character of 
strong interactions at high density is still not understood completely. At the core of protoneutron 
and neutron stars the matter density ranges from a few times of the density of normal nuclear 
matter to the value of one order higher than that when hyperons are expected to appear. The di- 
rect consequence of the extreme conditions inside protoneutron and neutron stars is the possibility 
of the appearance of different exotic forms of matter. Of special interest is the existence in these 



high density interiors the strangeness components 



characteristic mass-radius relation of the star [22j]. The appearance of the additional degrees of 



ike hyperons which may significantly change the 



freedom and their impact on protoneutron star structure and evolution have been the subject of 
extensive studies [a]. 

The properties of matter at extreme densities are of particular importance in determining forms of 
equations of state relevant to neutron stars and successively examining their global parameters P|. 
Theoretical description of hadronic systems should be performed with the use of quantum chromo- 
dynamics (QCD) as it is the fundamental theory of strong interactions. However, at the hadronic 
energy scale where the experimentally observed degrees of freedom are not quarks but hadrons the 
direct description of nuclei in terms of QCD becomes inadequate. In this paper an effective model 
)ased on chiral symmetry in its nonlinear realization has been introduced. The chiral SU(3) model 
which has been applied includes nonlinear scalar and vector interaction terms. This model 
offers the possibility of constructing a strangeness rich neutron star model and provides its detail 
description. 

Another alternative theoretical approach to the correct description of nuclear matter which has been 
formulated is quantum hadrodynamics (QHD) Q|- This theory gives quantitative description of 
the nuclear many body problem. QHD is a relativistic quantum field theory in which nuclear mat- 
ter description in terms of baryons and mesons is provided. The original model (QHD-I) contains 
nucleons interacting through the exchange of simulating medium range attraction a meson and u 
meson responsible for short range repulsion. The extension (QHD-II) of this theory includes 
also the isovector meson p. Nonlinear terms in the scalar and vector fields were added in order to 
get the correct value of the compressibility of nuclear matter and the proper density dependence in 
the vector self-energy. The variation of nucleon properties in nuclear medium is the key problem in 
nuclear physics. The mentioned above self-consistent relativistic field models involving coupling of 
baryons to scalar and vector mesons are succesfull in describing many properties of nuclear matter. 
However, there is a problem connected with the fact that the nucleon effective mass becomes very 
small at moderately high density The situation is even worse when hyperons are included. 
The onset of hyperon formation depends on the hyperon-nucleon and hyperon-hyperon interac- 
tions. Hyperons can be formed both in leptonic and baryonic processes. Several relevant strong 
interaction processes proceed and establish the hadron population in neutron star matter. When 
strange hadrons are taken into account uncertainties which are present in the description of nuclear 
matter are intensified due to the incompleteness of the available experimental data. The standard 
approach does not reproduce the strongly attractive hyperon-hyperon interaction seen in double 
A hypernuclei. In order to construct a proper model which would include hyperons the effects of 
hyperon-hyperon interactions have to be taken into account. These interactions are simulated via 
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(hidden) strange meson exchange: scalar meson /q(975) (a* meson) and vector meson ^(1020) ((f) 
meson) and influence the form of the equation of state and neutron stars properties. 
The solution of the presented models are gained with the mean field approximation in which meson 
fields are replaced by their expectation values. The parameters used are adjusted in the limit- 
ing density range around the saturation density po and in this density range they give very good 
description of finite nuclei. However, incorporation of this theory to higher density requires an 
extrapolation which in turn leads to some uncertainties and suffers from several shortcomings. The 
standard TM1 parameter set for high density range reveals an instability of neutron star matter 
which is connected with the appearance of neg ative nucleon effective mass due to the presence of 
hyperons. The Zimanyi-Moszkowski (ZM1) |23j model in which the Yukawa type interaction g a N& 
is replaced by the derivative one (g a N^r /MN)if)N^ u d u ifi exemplifies an alternative version of the 
Walecka model which improves the behaviour of the nucleon effective masses. It also influences 
the value of the incompressibility K of neutron star matter. The derivative coupling effectively 
introduces the density dependence of the scalar and vector coupling constants. 
The chosen models are very useful for describing properties of nuclear matter and finite nuclei. Its 
extrapolation to large charge asymmetry is of considerable interest in nuclear astrophysics and par- 
ticulary in constructin g pr otoneutron and neutron star models where extreme conditions of isospin 



onstructmg pre 



are realized [2, Liu LUj- LUj| ■ The models considered describe high isospin asymmetric matter and 
require extension by the inclusion of isovector-scalar meson ao (980) (the S meson) and nonlin- 
ear meson interaction terms. This extension affects the protoneutron stars chemical composition 
changing the proton fraction which in turn affects the properties of the star. 



I. PROPERTIES OF PROTONEUTRON STAR MATTER. 

The collapse of an iron core of a massive star leads to the formation of a core residue which 
is considered as an intermediate stage before the formation of a cold, compact neutron star. This 
intermediate stage which is called a protoneutron star can be described H, Q, Q, Q as a hot, 
neutrino opaque core, surrounded by a colder neutrino transparent outer envelope. The evolution 
of a nascent neutron star can be described by a series of separate phases starting from the moment 
when the star becomes gravitationally decoupled from the expanding ejecta. In this paper two 
evolutionary phases which can be characterized by the following assumptions: 

• the low entropy core s = 1 — 2 (in units of the Boltzmann's constant) with trapped neutrinos 
Y L = 0.4 
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I PROPERTIES OF PROTONEUTRON STAR MATTER. 



• the cold, deleptonized core {Yl = 0, s = 0). 

have been considered. These two distinct stages are separated by the period of deleptonization. 
During this epoch the neutrino fraction decreases from the nonzero initial value (Y v / 0) which 
is established by the requirement of the fixed total lepton number at Yj = 0.4, to the final one 
characterized by Y v = 0. Evolution of a protoneutron star which proceeds by neutrino emission 
causes that the star changes itself from a hot, bloated object to a cold, compact neutron star. 
The interior of this very early stage of a protoneutron star is an environment in which matter 
with the value of entropy of the order of 2 with trapped neutrinos produces a pressure to oppose 
gravitational collapse. The lepton composition of matter is specified by the fixed lepton number 
Yl = 0.4. Conditions that are indispensable for the unique determination of the equilibrium 
composition of a protoneutron star matter arise from the requirement of (3 equilibrium, charge 
neutrality and baryon and lepton number conservation. The later one is strictly connected with 
the assumption that the net neutrino fraction Y u / and therefore the neutrino chemical potential 
fi u 0. When the electron chemical potential // e reaches the value equal to the muon mass, muons 
start to appear. Equilibrium with respect to the reaction 

e~ + v e <r* fT + (1) 

is assured when /x M = fj, e — fj, Ue (setting fi u ^ = 0). The appearance of muons reduces the number of 
electrons and also affects the value of the proton fraction in matter. In the interior of protoneutron 
stars the density of matter can substantially exceed the normal nuclear matter density. In such a 
high density regime, it is possible that nucleon Fermi energies exceed the hyperon masses and thus 
additional hadronic states are expected to emerge. The higher the density the greater number of 
hadronic species are expected to appear. They can be formed both in leptonic and baryonic pro- 
cesses. The chemical equilibrium in stellar matter establishes relation between chemical potentials 
of protoneutron star matter components. In the case when neutrinos are trapped inside matter the 
requirement of charge neutrality and equilibrium under the week processes 



B 1 ^B2 + f + v f B 2 + f B x + v f 



(2) 



leads to the following relations 




(3) 



4 



A Scalar mesons. 



II THE MODEL. 



where 6j is the baryon number of particle i, (ft is its charge, / stands for leptons / = e, fi and 
(J>uf = Me- The nonzero neutrino chemical potential changes chemical potentials of the constituents 
of the system and also influences the onset points and abundance of all species inside the star. The 
increase of electron and muon concentration is a significant effect of neutrino trapping. The proton 
fraction Y p also takes higher value in order to preserve charge neutrality. In the case of strangeness 
rich matter the appearance of charged hyperons permits the lower electron and muon content and 
thus the charge neutrality tends to be guaranteed with the reduced lepton contribution. After 
deleptonization the neutrino chemical potential reduces to zero and chemical equilibrium inside 
matter differs from the one presented above. The relations between chemical potentials of different 
constituents of the system is now obtained setting \x v = 0. This stage is followed by an overall 
cooling stage during which the entropy in the star decreases. 



II. THE MODEL. 



The most general form of the Lagrangian function C of the whole system can be shown as a sum 
of separate parts which represent baryon together with baryon-meson interactions, scalar, vector 
and lepton terms, respectively 

C = C B + C s + Cv + C L . (4) 



A. Scalar mesons. 



The construction of the U(3)l x U(3)r invariants [3|,|2l| allows us to include different forms of 
meson-meson interactions. The following terms have been introduced 

X = Tr($ t $) ! Y = Tr[($ f $) 2 ], Z = det$> + det&. (5) 

where the first two are U(3)l x U(3)r invariant, whereas the Z term explicitly breaks the U(1)a 
symmetry. The presented above invariants are constructed from the 3x3 matrix field $ which 
enables the collective representation of the spin fields 

$ = ^Tafa = ^=Aa(oo + in a ), (a = 0, . . . , 8) (6) 



where T a = X a are the generators of U(3) and A a are the Gell-Mann matrices with Ao = y |l 
included for the singlet. The A matrices are normalized by 

Tr(A a A a ) = 25 ab . (7) 
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The <7 a and 7r a fields are members of the scalar £ = -^^a&a (J P = + ) and pseudoscalar II 
;^A a 7r a (J p = 0~) nonets respectively. Presenting them as matrices one can obtain 



V2 L 



O"0 



75 a ° + 76 a8 + 73 ao 



1 _L 1 ~ 



(8) 



n 



(9) 



The scalar part of the Lagrangian function £5 is a sum of the U(3)l x J7(3)# symmetric term and 
the explicit symmetry breaking one 

Cs(<S>)=£s,sy m ($) + £ S B- (10) 
The £s,sym(&) part in turn contains the kinetic and potential terms and can be expressed as: 
£s, sym ($) = ^Tr(D^D^) - [/($) (11) 
with the potential function U(&) defined in the following way 

[/($) = m 2 $t$ + ai [Tr(&<S>)} 2 + a 2 Tr($ t $) 2 (12) 

where m is the tree level mass of the fields in the absence of symmetry breaking term, a\ and 0:2 
are coupling constants. The explicit symmetry breaking term Csb has the following form: 

C SB = c[Det{<&) + Det(&)] + Tr[H($ + (13) 

The first term which breaks the U{1)a symmetry of the Lagrangian explicitly gives the mass 
to the pseudoscalar singlet. In the second term H denotes the 3x3 matrix with H = 
\j ^[2diag(m\^,m\^,2m\j K - m 2 /^) where m n = 139 MeV, m K = 498 MeV. 
The spontaneous chiral symmetry breaking is triggered by a non-vanishing expectation value a 
corresponding to the location of the minimum of the potential U(&). After the symmetry breaking 
the $> field acquires a vacuum expectation value. Owing to parity conservation in infinite neutron 
star matter the pseudoscalar fields ir a cannot assume a non-vanishing vacuum expectation value 
for <£. Shifting the scalar fields by their vacuum expectation values and substituting them to the 
Lagrangian function (jTU)l allows us to determine the potential function £/(<£— < >). 



A Scalar mesons. 



II THE MODEL. 



Ttte Potential 




Figure 1: The form of the potential function U(ip, (). 

The scalar fields a a in the basis of U(3) generators are not mass eigenstates and the obtained 
mass matrix is not diagonal. Since the mass matrix is real and symmetric, there exists a real 
orthogonal matrix O which diagonalizes M.. On the assumption that cp = (1/V6a$ + l/v^o) 
and C = + l/V^o), the physical scalar fields are obtained as a result of the following 

procedure 



cos-d siwQ 
—siwQ cos-d 





M 






)i 




N 






M 







(14) 



which yields the diagonalization and shifting the fields to the minimum of the potential U(ipX)- 
The form of the potential function U((p,Q is depicted in Fig^ For the remaining a{] state the 
following substitution has been established a[| = 5. The field a describes a broad resonance {ra a ~ 
470 MeV) connected with the exchange of a correlated pair of pions. 

B. Baryon-meson interaction. 



Baryon fields that enter the model are grouped into a 3 x 3 traceless hermitian matrix 



B 



\ 



(15) 



For the baryon kinetic terms to preserve chiral invariance the local, covariant derivatives have to 



be used 
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24] 



D^B = d fl + i[T ll ,B] 



(16) 
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where 

T^ = - l -Wd^u + ud^). (17) 
and 

i 

u = exp[— 7r a A a 7 5 ]. (18) 

Thus the pseudo-scalar mesons are given as parameters of the symmetry transformation. 
The general form of the baryon(-B) meson (W) interaction terms can be written as follows 

Cbw = -V2gf {a w \BOBW] F + (1 - a w )\BOBW] D ) - gf ' -=Tr(BOB)TrW (19) 

V3 

where 

[BOBW] F = Tr(BOWB - BOBW) (20) 

[BOBW] D = Tr(BOWB + BOBW) - -Tr(BOB)TrW (21) 

The relation ffTTTf) it is a mixture of the F-type (antisymmetric) and D-type (symmetric) couplings, 
aw denotes the F/F + D ratio. The differences for the baryon-scalar and baryon- vector meson 
interactions are connected with difference in Lorentz space. For scalar mesons W = $ and O = 1 
The masses of the whole baryon multiplet are generated spontaneously by the vacuum expectation 
value of the nonstrange and strange meson condensates 

<M >=diag(^=,^=,() (22) 

The assumption that ay = 1 has be made. Setting also g\ = V6g^ the model in which nucleon 
mass does not depend on the strange condensate is obtained. The value of the nucleon mass can 
be obtained by adjusting only one coupling constants. The correct values of the remaining baryons 
it is necessary to introduce an explicit symmetry breaking term. 

For the baryon-vector meson interactions W = and O = 7 M The universality principle and 
the vector meson dominance model points to the conclusion that the D-type coupling should be 
small. Setting as for the case of scalar mesons ay = 1 and g\ = v6ff8 , the latter assumption 
corresponds to the case when strange vector field does not couple to nucleon, the following form of 
the baryon-vector meson interaction Lagrangian can be written 

C BV = -V2g%([B^BV»] F + Tr(B^B)TrS^). (23) 
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with the singlet vector meson state and being the matrix containing the octet of vector meson 
fields 

\ 



( J_T/8,J_ 



V 



K 



J_T/8 _ J_ 

K*° 



if* 



K*° 



%/6 V / 



(24) 



The relations for the vector meson couplings reflect quark counting rules. 

The physical meson states u> and 4> are mixed states. They stem from the singlet and octet 
vector meson states and are given by the relation 



lo = cos 6 — sin ( 



(25) 



V°sm9 + S^cosd. 



Under the condition that the vector meson (ft is nearly a pure ss state as it decays mainly to kaons 
the mixing with the mixing angle tan 9 = l/y/2 is called ideal. 



C. Vector mesons. 



The general vector meson Lagrange function can be presented as a sum of the kinetic energy 
term, mass term and higher order self-interaction terms: 

C v = --TriV^Vn + ^m 2 v Tr(V^) - c[Tr{V^f] ~ dTr[{V^) 2 }. (26) 

with the field tensor V^ u = d^V v — d v V^. This form of the mass term implies a mass degeneracy 
for the meson nonet. 



D. Leptons. 

Leptons are gathered into SUl(2) dublets and singlets with the flavor / (/ = e,/x) 
/ 

, e fR , v fR 

L 



Assuming that neutrinos Vf have masses m/ the lepton Lagrange function has the form 
~ E /,/' M f:f ,(v fL v fR + h.c.) - E/,/' Mfj,(V fL v f , R + h.c.) 
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The physical massive neutrinos are mixture of neutrinos with different flavors 

VfL = Ei U LJi u iL 
u fR = J2i Unjoin 

III. RELATIVISTIC MEAN FIELD EQUATIONS. 

To investigate the properties of infinite nuclear matter, the mean field approximation has been 
adopted. The symmetries of infinite nuclear matter simplify the model to a great extent. The 
translational and rotational invariance claimed that the mean fields of all the vector fields vanish. 
Only the time-like components of the neutral vector mesons have a non-vanishing expectation 
value. Owing to parity conservation, the vacuum expectation value of pseudoscalar fields vanish 
(< Ti a >= 0). Meson fields have been separated into classical mean field values and quantum 
fluctuations, which are not included in the ground state. Thus, for the ground state of homogeneous 
infinite nuclear matter quantum fields operators are replaced by their classical expectation values. 



a = a +sq a* = a* + s* 5 a = 5° + d S 3a 

4>ll = 4>v + /o^o ^ = uJ M + woS^ b^=V^ + r <^ <5 3a 

When strange hadrons are taken into account uncertainties which are present in the description of 
nuclear matter are intensified due to the incompleteness of the available experimental data. The 
standard approach does not reproduce the strongly attractive hyperon-hyperon interaction seen in 
double A hypernuclei. In order to construct a proper model which do include hyperons the effects of 
hyperon-hyperon interactions have to be taken into account. These interactions are simulated via 
(hidden) strange meson exchange: scalar meson /o(975) (a* meson) and vector meson </>(1020) ((f) 
meson) and influence the form of the equation of state and neutron stars properties. The presence of 
hyperons demands additional coupling constants which have been fitted to hypernuclear properties. 
The Lagrangian function for the system is a sum of a baryonic part, including the full octet of 
baryons, baryon-meson interaction terms, a mesonic part and a free leptonic one. The Lagrangian 
density of interacting byarons jCbm can be explicitly written out in the following form: 

CbM = i i'Bl^D^B - Y ^B M eff,B^B (27) 
B B 
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where the spinor ^ T = V'Ai "0£> V'h) is composed of the following isomultiplets 



N 



* 5 



_D M is the covariant derivative given by 



(28) 



The meson part of the Lagrangian function includes nonlinear terms that describe additional inter- 
actions between mesons. 



Cm = ^ad^a+hs^d^ + U^d^o* - U effjS (a,6,a*) 



1 



1 



(29) 



The field tensors Q^,^^ and i?" are defined as 



R 1u — 9fj.pl - d u p® + g p €-abcp\p C v 



(30) 
(31) 



All nonlinear meson interaction terms are collected in the vector and scalar potential functions 

U e ff tV (u,p,(j)) and U e ff tS (a,S,a*). 

The form of the scalar potential is portrayed by 



U eff , s (a,5,a*) = \ir?y + ±g 2 a 3 + \g^ + \m 2 5 {5 a 5 a ) 
+\gsA{S a 5 a f + g a5 2<r5 a 6 a + \g u2 52'J 2 5 a 5 a 



1 2 *2 1 



*3 



1 



*4 



*2 



2 m a* cr +^9a*3CT + -g a *AV + g<jo*2VO- 
+SW*3°" cr * 3 + 9a2a*C 2 ^* + 5(72(7* 2 cV* 2 + g r7 3cr* a ' 3 

+g82a*5 2 <J* + gs2a*2S 2 a* 2 + g a 82a*o-5 2 a* , 



(32) 



whereas the effective vector potential has the following form 

U eff , v (co, p, cj>) = \c 3 (K^) 2 + {plp^) 2 - 3(w/ + p> a ")(<^) 



^ 2 



(33) 
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Effective Potential U e(( (o, 0, 0) 




Figure 2: The scalar potential function calculated for different parameter sets. 



Table I: Parameter sets for chosen equations of state. 





TM1 


GM3 


SU3 


m a (MeV) 


511.2 


450.0 


477.6 


m s (MeV) 


980.0 


980.0 


1029 


92 


7.2325 fmT 3 


15.286 fm' 3 


16.196 fm- 3 


93 


0.6183 


-6.4547 


15.669 


9<tN = 9(7 


10.029 


7.186 


7.923 


9lu 


12.61 


8.702 


9.225 


9p 


9.264 


8.542 


5.742 


9aS2 








32.639 fm-' 3 


9a252 








29.511 


95A 








19.781 



The potent taction V« lM i S presented in F lg0 Scat, potential functions generated for 
TM1 and GM3 (Glendenning and Moszkowski [9]) parameter sets are also included in this figure. 

The scalar a coupling constants for hyperons are chosen to reproduce hyperon potentials in 
saturated nuclear matter: 



30MeV, U£ (p ) = +30M e y, (p ) = -\8MeV 



A' / 



(34) 



Upon a recent analysis of S~ atomic data there has been indication of the existence of a repulsive 
isoscalar potential in the interior of nuclei. This corresponds to a resent search for £ hypernuclei 
which show the lack of bound-state or continuum peaks. The only found E bound state is f^He. 
The binding results from the strong isovector component of the X nuclear interaction. Thus, in 
the considered model X hyperons wil not appear in the neutron star interior. The experimental 
data concerning hyperon-hyperon ^| interaction are extremely scarce. The observed double A 
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Table II: The nonlinear SU3 parameter set. 



m CT (MeV) 


m 5 (MeV) 


92 


93 


9a 


477.6 


1029 


16.196 fm- 3 


15.669 


7.923 


9u 


9 P 


95A 


9a52 


9a282 


9.225 


6.754 


19.78 


31.69 fm- 1 


59.02 




982a* 


982a*2 


9aS2<T* 


9 17*4 


49.88/m- 1 


-9.51 


1.90 


21.39 


34.07 






9o2a* 


50-20-* 2 


9cr3cr* 


25.27 


12.01 


-9.51 


1.90 


4.88 



hypernuclear events in emulsion require strong AA attractive interaction. An analysis of events 
which can be interpreted as the creation of AA hypernuclei allows us to determine the well depths 
of hyperon in hyperon matter. The a* coupling constants to hyperons, in accordance with the 
one-boson exchange model D of the Nijmegen group and the measured strong AA interaction, are 
fixed by the condition: 

- « -£/{ S) « -2[/i A) « -2f/| A) « 40Af e y. (35) 

The appropriate parameter set is denoted as SU3 constrained not only by the value of physical 
scalar meson masses but also by the properties of nuclear matter at saturation. The obtained 
results are collected in Table 1 and 2. 

The effective interaction is introduced through Klain-Gordon equations for the meson fields with 
baryon densities as source terms. These equations are coupled to the Dirac equations for baryons. 
All the equations have to be solved self-consistently. The field equations derived from the Lagrange 
function at the mean field level are the following: 

(ml + g a 252dl + ga2a*S* )s + (g 2 + ga3a*S* )sl + 5 3 S(j ( 36 ) 

= J29vBM* ff S(M eff , B ) 

B 

3 

mlw + c 3 Wo - -c 3 w fo = g^BUB (37) 

1 B 

in 2 p r + c 3 ro - -c 3 r /o = g P BhBn B (38) 

1 B 



mjdl + 2(g aS 2So + 7; 9*262 si + -rg54 + 552o*So + g52**2sf + ga62a*S* S )d (39) 
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^9SBhBS{M e ff )B ) 



B 



+ g52a*2dl)s* (40) 

+ (9a2a*sl + 9a3cT*sl+ g82a*dl) = ^ g a * B Mjf fS(M e ff >B ) 

B 

3 

m lfo ~ ^c 3 (wl + ro)/ - c 3 /o = 94>sn B - (41) 

1 B 

The function S(M e ff )B ) is expressed with the use of the integral 

S( M«„, B) = ^f —^L-tfo-fc, (42) 

where </# and I3B are the spin and isospin projection of baryon B, n B which denote the baryon 
number density is given as 
2J B + 1 



n B 



POO 

/ k 2 dk(f B -f § ) (43) 
Jo 



2vr 2 

The functions f B and f B are the Fermi-Dirac distribution for particles and anti-particles respectively 

= 1 + e E[(k,Meff,B)T^B]/k B T- ( 44 ) 

The Dirac equation for baryons that is obtained from the Lagrangian function has the following 
form: 

{i-fdp - M B>eff - guB^wo ~ g^Bl°fo - ^g P Bl°T 3 r )i; = (45) 

from this equation it is noticeable that the baryon develop the effective mass M Be ff generated by 
the baryon and scalar field interactions and defined as: 

M Bjf ,ff = M B - (g aB s Q + g a * B s* Q + h B gs B do) (46) 
and the effective chemical potential ix e ff is given by 

t*B,eff = HB- (gwBl°^0 ~ g<t>Bl Q h - ^9pBl°T 3 r ) (47) 

Numerical solutions of the equation lj4T?jl is presented in FigEl The main effect of the inclusion 
of S meson itself becomes evident when studying properties of neutron star matter especially 
baryon mass splitting and the form of the equation of state. FigEI depicts the effective baryon 
masses as a function of baryon number density n B . There is a noticeable mass splitting of 
baryons for each isomultiplet in medium which depends on the considered baryonic masses. The 
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1000 r i | i | i | i | i | i | i | i | i | i | i | i 




Figure 3: The effective nucleon masses obtained for different parameter sets. 



Table III: Meson-hyperon coupling constants for the ZM1 model. 
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bigger is the mass, the smaller the mass differences are. Depending on the sign of the third 
component of particular baryon isospin the S meson interaction increases the proton and H° 
effective masses and decreases masses of neutron, S~ and E~. Contrary to this situation, when <5 
meson is not included, the baryon mass for a given isomultiplet remains degenerated. Through- 
out the effective baryonic masses 5 meson alters baryon chemical potentials what is reflected in 
characteristic modification of the appearance, abundance and distributions of the individual flavors. 



As it was stated earlier neutron star matter possesses highly asymmetric character caused by the 
presence of small amounts of protons and electrons. The introduction of the asymmetry parameter 
f a which describes the relative neutron excess defined as 

fa = (48) 

n N 

allows to study the symmetry properties of the system. Similarly to the asymmetry parameter /„ 
a parameter f s which specify the strangeness content in the system and is strictly connected with 
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the appearance of particular hyperon species in the model has been introduced. 

f s = nA + nS + 2n5 . (49) 
n\ + ™£ + n~ + nN 

IV. THE DERIVATIVE COUPLING MODEL. 



The description of the nuclear matter properties with the use of the standard TM1 parameter 
set reveals a shortcoming which is connected with the appearance of negative nucleon effective 
masses for densities characteristic for hyperon stars. The Zimanyi-Moszkowski (ZM1) j^J model 
in which the Yukawa type interaction g a N<J is replaced by the derivative one (g a N& /Mn^nIvQ 1 'ifi 
exemplifies an alternative version of the Walecka model and improves the behavior of the nucleon 
effective masses. 

In this model in the term which represents the Lagrangian density of interacting baryons is given 
as follows 

r S^n , 9aBa + 9a* B<?* + hB9SBT a 5 a - . 

B B B 



Rescaling the baryon field in a way proposed by Zimanyi and Moszkowski [2j| the modified Lagrange 
function for interacting baryons is obtained 

r v-7 ■ , srft , 9uB(y + 9a*B(y* + hB95BT a 5 a \~ 1 - 

C-B = -l^VBlTDnip -Z^( 1 + m ) ^B M B^B (51) 

B B B ' 

Expanding the expression (1 + g a N(T /M^)^ 1 in terms of g s Ncr/M]^ (index N denotes nucleons) up 
to first order in a allows one to reproduce the baryonic part of the Lagrangian of the Walecka model 
with a Yukawa N — a interaction. The effective baryon mass in this case is given by the following 
formula 

M R ff = 1 (52) 

"> e}} l+(gaB<T + go*B<T*+hBg5 B T a d a )/M B v ' 

The parameters employed in the ZM1 model are collected in Tables 3 and 4 jjj. In this case 
the value of the parameter g p N has to be redefined in comparison with the standard TM1 value. 
The parameters g p N and gsN are adjusted to obtain the symmetry energy coefficient a sym (no) at 
saturation equal 32.4 MeV which is in good agreement with the empirical value of about 30 ± 4 
MeV. 
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Table IV: Meson-nucleon coupling constants for the ZM1 model. 
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V. RESULTS AND CONCLUSIONS. 

On obtaining the form of the equations of state of protoneutron star matter the mass-radius 
relation and composition of the star can be specified. Two distinctive parameter sets describing 
strangeness rich matter have been used. The first one denoted as SU3 has been constructed on 
the basis of chiral symmetry whereas the latter one (ZM1) has been obtained with the use of the 
derivative coupling model. The considered theories have been extended by the inclusion of S meson 
and nonlinear vector meson interaction terms. The inclusion of 5 meson seems to be indispensable 
for the complete description of asymmetric neutron star matter. The main effect of the presence of 
5 meson becomes evident when studying properties of matter especially baryon mass splitting and 
the form of the equation of state. Throughout the effective baryon masses 5 meson alters baryon 
chemical potentials what realizes in characteristic modification of the appearance, abundance and 
distribution of the individual flavors. The assumptions underlying the calculations performed in 
this paper are connected with the choice of the repulsive nucleon-hyperon S interaction. 
The forms of the selected equations of state are shown in Fig0J In agreement with the generally 
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Figure 4: The pressure energy density relations for the selected parameter sets. 
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Figure 5: The equilibrium composition of protoneutron star matter as a function of the baryon number 
density tlb- Upper left panel depicts the composition obtained on the basis of the ZM1 model for s = 2 and 
Yi = 0.4. Upper right has been calculated for s = and Yi = 0. Results obtained for SU3 model are presented 
in lower panels s — 2 and Yi — 0.4 (left) and s — and Yi — (right). 

accepted scheme of protoneutron star evolution results for two different cases has been presented. 
The first one corresponds to the era of neutrino trapping (Yl = 0.4 and s = 2) and the second one 
fulfills the conditions of cold, deleptonized matter (Yl = and s = 0). This figure presents also 
the influence of both entropy and neutrino trapping. For these two cases the SU3 models lead to 
significantly stiffer equations of state. 

In the four successive panels of FigEl the analysis of the relative concentrations of particles as 
functions of baryon number density rig are presented. These results indicate that the first strange 
baryon that emerges is the A it is followed by E~ and H°. For both parameter sets the sequence of 
appearance of hyperons is the same however, in case of ZM1 parameter set, shifted towards higher 
densities. Due to the repulsive potential of S hyperons their onset points are possible at very high 
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model 


M{M S ) 


R(km) 


p c (10^ g / cm 6 ) 


T c (MeV) 


P c (MeV/fm 3 ) 


SU3, S = 2, Y L = 0.4 


1.88 


12.3 


2.0 


49.68 


268.1 


SU3, T = 


1.49 


10.52 


2.43 


0.0 


189.9 


ZMl, S = 2, Y L = 0.4 


1.62 


11.41 


2.44 


34.88 


300.4 


ZMl, T = 


1.46 


11.2 


1.98 


0.0 


166.0 



Table V: Properties of protoneutron maximum mass configurations for SU3 and ZMl parameter sets. 



densities which are not relevant for neutron stars. The appearance of H~ hyperons through the 
condition of charge neutrality affects the electron and muon fractions and causes a drop in their 
contents. Therefore the appearance of charged hyperons permits the lowering of lepton contents 
and charge neutrality tends to be guaranteed without lepton contribution. Larger effective baryon 
masses cause the shift of given hyperon onset point especially for the charged ones in the direction 
of higher densities. 

Any realistic calculation of the properties of neutron stars is based upon the general relativistic 
equation for hydrostatic equilibrium (the Oppenheimer-Volkoff equation j^). 
On specifying the equation of state the solution of this equation can be found and this allows us 
to determine global neutron stars parameters Q|. The structure of spherically symmetric neutron 
star is determined with the use of this equation which is integrated starting from p = p c at r = 
to the surface at r = R where p = 0. In this way a radius R and a mass M = m{R) for a given 
central density can be found . As a result, the mass-radius relation can be drown. These relations 
for the chosen forms of the equations of state are presented in FigEJ This particular sequence of 
figures gives a detailed inside into what happens with lowering number of neutrinos (Y v — > 0) and 
decreasing value of entropy. In the left panel of FigEJthe mass-radius relations for the protoneutron 
star calculated with the use of ZMl parameter set is plotted. The right panel depicts the mass- 
radius relations for SU3 model. In general neutrino trapping increases the value of the maximum 
mass. The SU3 parameter set gives higher value of the maximum mass then the ZMl parameter set. 
Dots which are connected by straight lines represent the evolution of configurations characterized 
by the same baryon number. In both figures there are configurations which due to deleptonization 
go to the unstable branch of the neutron star mass-radius relation. The reduction in mass is 
much bigger in case of the SU3 parameter set. The stellar parameters obtained for the maximum 
mass configuration for different models are presented in Table The solutions of the structure 
equations also allows us to carry out an analysis of the onset point, abundance and distributions 
of the individual baryon and lepton species as functions of the star radius. The comparison of 
results obtained for the two cases presented above (ZMl and SU(3)) have been made on the basis 
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Figure 6: Deleptonization of a protoneutron star. The left panel presents results obtained for ZM1 parameter 
set and the right panel for the SU3 model. 

of the assumption that the repulsive £ interaction shifts the onset point of £ hyperons to very 
high densities. Consequently, they do not appear in the neutron star interior calculated in these 
models. The maximum mass configurations has been considered. In case of trapped neutrinos the 
ZM1 model leads to more compact strangeness rich core. This very compact hyperon core which 
emerges in the interior of the maximum mass configuration consists of H~ and A hyperons as shown 
in the upper left panel of FigQ The presence of negatively charged hyperons reduces the content 
of negatively charged leptons. For the SU3 model the strangeness rich core is more extended (the 
bottom left of FigEJ but in case of neutrino trapped matter it contains only A hyperon. The 
absence of negatively charged hyperons leaves the electron and muon content unchanged. After 
deleptonization the strangeness rich core in case of the SU3 model remains more extended than in 
the case of ZM1 parameter set but it contains E°, H~ and A hyperons. 

The relative baryon composition in this model can be also analyzed through the density dependence 
of the asymmetry parameter f a and the strangeness contents f s . FigJHl presents both parameters as 
functions of the star radius R. As it was shown in Figd trapped neutrinos have an influence on the 
charged hyperon onset points. The appearance of additional negatively charged particles has the 
consequence on the proton content in the system and this in turn on the asymmetry parameter f a . 
Fig. El depict the relative fractions of protons, electrons and muons for the ZM1 and SU3 parameter 
sets. Fig. [HI compares the asymmetry parameter f a for ZM1 and the SU3 parameter sets. It is 
noticeable that SU3 parameters give less asymmetric protoneutron and neutron star but with the 
higher value of the strangeness content parameter f s . 
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Figure 7: The equilibrium composition for the maximum mass configuration as functions of the stellar radius 
R. Upper left panel depicts the composition obtained on the basis of the ZM1 model for s — 2 and Y\ = 0.4. 
Upper right has been calculated for s = 0. Results obtained for SU3 model are presented in lower panels 
s = 2 and Yi = 0.4 (left) and s = (right). 

Fig. depicts the temperature profiles plotted for the presented above models. In case of ZM1 
protoneutron star model the temperature in interior of maximum mass configuration takes higher 
value than in case of SU3 model. The appearance of the strangeness-rich core on both plots is 
marked by the local inversion of temperature. 
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panel) as functions of the stellar radius R . 
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